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ABSTRACT 

We present period measurements of a large sample of field stars in the solar neigh- 
bourhood, observed by CoRoT in two different directions of the Galaxy. The presence 
of a period was detected using the Scargle Lomb Normalized Periodogram technique 
and the autocorrelation analysis. The assessment of the results has been performed 
through a consistency verification supported by the folded light curve analysis. The 
data analysis procedure has discarded a non-negligible fraction of light curves due to 
instrumental artifacts, however it has allowed us to identify pulsators and binaries 
among a large number of field stars. We measure a wide range of periods, from 0.25 
to 100 days, most of which are rotation periods. 

The final catalogue includes 1978 periods, with 1727 of them identified as ro- 
tational periods, 169 are classified as pulsations and 82 as orbital periods of binary 
systems. Our sample suffers from selection biases not easily corrected for, thus we do 
not use the distribution of rotation periods to derive the age distribution of the main- 
sequence population. Nevertheless, using rotation as a proxy for age, we can identify a 
sample of young stars 600 Myr), that will constitute a valuable sample, supported 
by further spectroscopic observations, to study the recent star formation history in 
the solar neighborhood. 

Key words: stars: solar-type - stars: rotation - solar neighbourhood. 



1 INTRODUCTION 

CoRoT is a pioneering space mission that provided a large 
scale survey of photometric variability of stars with a very 
high photometric precision and long time coverage. CoRoT 
observed a large sample of stars consisting primarily of 
main sequence stars in the solar neighbourhood. CoRoT 
primary goals are stellar seismology and extrasolar planets 
search, nevertheless, stellar light curves are themselves a rich 
source of astrophysical information, thanks to their luminos- 
ity modulation. There are several reasons for which a star 
can vary its luminosity, and these may be linked to starspot 
modulation due to stellar magnetic activity, to low and high 
amplitude pulsations, to the existence of close companions 
(binary systems), to the periodic crossing of exoplanets in 
front of the disk of the star, and so on, including a complex 
mixing of several of the formerly mentioned phenomena. Of 
particular interest is the study of the periodic modulation 
of starlight produced by non-uniformities on the surface of 
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a main-sequence star due to manifestation of stellar activ- 
ity (e.g., spots and plages) which are used to determine the 
stellar rotation period. 

Stars contract during the pre-main sequence phase, but 
the mechanisms governing the angular momentum evolution 
during contraction are still not well understood (depending 
basically on the star-disk interaction). They reach the zero 
age main sequence (ZAMS) with a maximum velocity, and 
stars of the same spectr al type are observed to have a spread 
of rotational velocities (iB ouvier et al.|[l997l ). 

During the main sequence (MS) phase rotation de- 
creases with age (IWilsonf [196^ : lKraft|]l967l) . due to the 
loss of angular m omentum through magnetized stellar wind 
(|Schatzmanlll96 j ). and all the F, G and K-type stars tend to 
converge to a single, color-dependent, rotation-age relation 
for stars older than ~ lO^yr corresponding approximately to 
the age of the HyadesQ Therefore, when a star reaches this 
stage, stellar rotation can be used in principle to estimate 
stellar ages. The first author that determined a relat ion be- 
tween rotation and age has been ISkumanichI (|l972l ) which 
obtained Prot oc where t is the stellar age. More 



^ The described scenario is valid for single stars, while tidally 
locked binary systems maintain fast rotation even at old age. 
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recently more complex relations, called gyrochronology re- 
lations, that ta ke into account also the stellar mass, have 
been proposed (|Barnesll200^ . l2007l : iMamaiek fc HiUenbrandl 
|2008| ). As consequence single stars with rotation faster than 
the convergence value cannot be precisely dated but by them 
we can put an upper limit of 10^ yr to their age, while slower 
stars of a given mass can be dated just from their rotation. 

Using CoRoT photometry we are able to reveal lumi- 
nosity variation, with a precision down to 0.1 mmag per 
hour (magnitude between 11 and 16), during continuous ob- 
servations (up to more than 150 days), allowing to measure 
photometric periods also in relatively quiet stars (for com- 
parison, the luminosity variations of the Sun range between 
~0.3 mmag and ~0. 07 mmag at maximum and minimum 
activity, respectivelv. I Aigrain et al. 2004). 

The search for periodic modulation changes in the light 
curves of dwarf stars observed by CoRoT will allow us to 
measure stellar rotation periods in a large sample of field 
stars, thus enabling an analysis of the period as function of 
spe ctral type (or stell ar mass). 

iBasri et al.l (|201ll ') recently performed a similar analysis, 
discussing the photometric variability characteristics in Ke- 
pler target stars, with a higher photometric precision. They 
preliminarly conclude that a large number of dwarf stars do 
vary because of rotation and starspots, and that the stars 
whose variability is periodic but not due to spots (pulsators, 
eclipsers) are far less numerous and can be recognized with 
good efficiency, but warning that further analysis is required 
to distinguish between various sources of variability in more 
ambiguous cases. 

The immediate goal of the present study is the compila- 
tion of a catalogue of periods of main-sequence stars in the 
CoRoT fields. This catalogue will be analyzed in order to 
study the rotation properties of solar-like stars in the solar 
neighbourhood. In particular we are interested in identify- 
ing the fast rotator stars that, as explained above, are the 
younger population in the solar neighborhood. This sample 
will be very valuable for follow-up observations aiming at 
reconstructing the star formation history close to the Sun. 
The present paper is organized as follows: in Sects. [2] and 
[S] we describe the CoRoT observations and the reduction 
of the light curves. In Sect.|4]we describe the methods used 
to derive periods, in Sect. [S] we show the identification of 
rotational periods and the compiled catalogue. In Sects. [6] 
and[7l we discuss and summarize our findings. 



2 COROT OBSERVATIONS 

CoRoT (convect ion Rotation and planetary Transits, 
iBaglin et al.|[200^ ) is a space mission dedicated to stellar 
seismology and to the search of extrasolar planets. Besides 
the two main scientific programs, several other "additional 
programs" are made possible by its high photometric per- 
formances and its observing runs covering up to five months 
without interruption. The data presented here have been ob- 
tained in the context of the program "An unbiased study of 
rotation and stochastic variability and fiaring in all CoRoT 
target stars" (P.I. F. Favata). 

The CoRoT instrument consists of a telescope with an 
aperture of 27 cm. The field of view is a square of 2.8° x 2.8° 
covered by four 2048 x 2048 pixel EEV CCDs, two special- 



ized for the seismology mission, the other two for the exo- 
planet mission. CoRoT is restricted to point within two 10 
degree radius circles, which were selected at the intersections 
of the ecliptic and Galactic planes, one towards the Galactic 
anti-center (RA2ooo=06''50'"25^ Dec2ooo = -01°42'00") and 
the other towards the Galactic center (RA2ooo = 19''23'"34'', 
Dec2ooo=00°27'36"). 

The exoplanet CCDs are read every 32 s, and photom- 
etry for up to 6000 stars per CCD (12 000 per run) in the 
range 11.5 < R < 16 is performed on-board. Light curves 
with a sampling of 512 s are recorded (32 s sampling is 
available for up to 500 stars per CCD). For almost all stars 
brighter than i?=15, the fiux is separated into broad-band 
red, green and blue channels, containing ~ 40, 30 and 30% of 
the flux, respectively (the separation between the red, green 
and blue bands is adapted for each target stars). 

The raw data, as downloaded from the satellite (la- 
belled NO data) are further processed on the ground, with 
the CoRoT pipeline which corrects several effects. Amongst 
these are the electronic offset, gain, electromagnetic inter- 
ference, and outliers. The pipeline includes background sub- 
traction and partial jitter correction for the three colour 
channels, and also flags data points collected during the 
South Atlantic Anomaly (SAA) crossing. Th e results o f 
these corrections are labelled Nl data (Sam adi et al.ll2007l '). 
Further corrections of Nl data, including merging of 512 
s and 32 s samplings, calculation of heliocentric date, hot 
pixels detection (d ata are flagged) lead to the N2 data 
(ISamadi et al ]|2007'). The N2 data are the science grade data 
on which all scientific analysis is performed. The light curves 
are normally sampled at a rate of 512 s or oversampled at 32 
s; in several cases (in particular for chromatic light curves) 
the oversampling mode was triggered after the start of the 
run. 

In the present paper we analyzed Corot's first long run 
observations of the exoplanet field. The first long run ob- 
servations (LRcOl) started with a pointing direction close 
to the center of the Galaxy and lasted from May 16 to Oc- 
tober 15, 2007 (155 days). The first long run observations 
with a pointing direction close to the anti-center (LRaOl) of 
the Galaxy lasted from October 24, 2007 to March 3, 2008 
(130 days). The light curves obtained are nearly continuous 
with only a small number of gaps resulting from the periodic 
crossing of the SAA. 

In Figures[T]and [2]we show two examples of N2 data, for 
quiet and active stars with a binning of two hours. The long 
period trend which is evident in the two examples shown 
(and in almost all raw light curves) has been removed as 
described later in Sec. [3] 

The parent CoRoT sample was selected as a magnitude- 
limited sample with objects affected by crowding removed. 
The targets observed by CoRoT were selected using the in- 
form at ion g at hered in the EXODA T database (|Deleuil et al.l 
l2006l . l2009l : iMeunier et al] |2007| ). buih with dedicated 
ground based photometric observations in the visible. Due to 
the requirements of the exoplanet search it was important 
to select the dwarf population among all potential target 
stars, in order to increase the chance of detecting plane- 
tary transits. Therefore the resulting parent population is 
biased toward dwarf stars e ven if a substa ntial fraction of 
giants may still be present (| Aigrain et al.ir2009 ). An addi- 
tional, intrinsic bias present in magnitude-limited samples 
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Figure 1. Example of a raw light curve for a quiet star from the first CoRoT run towards the Galactic anti-center, with a binning of 
2h, showing an instrumental long term variation in the left panel, detrendcd and normalized in the right panel. Time is in days from 
2000 January 1. 
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Figure 2. Same as Fig.[T] but for an active star. 



is the increased fraction of older (more luminous and there- 
fore more evolved) stars and binaries (again more luminous 
than single stars of the same age). The effects of these biases 
on our results will be discussed in Sec. [G] 



3 SAMPLE AND DATA REDUCTION 

We initially selected from the CoRoT database (Archive 
of the COROT Data Center at Space Astrophysics Insti- 
tute, Orsay-France, http://idoc-corot.ias.u-psud.fr/), obser- 
vations of dwarf stars (labelled as luminosity class V) with 
colour index B — V between 0.44 and 1.4 performed during 
the LRaOl and LRcOl runs of the exoplanet field (the lumi- 
nosity class and colour information are from both the EXO- 
DAT/COROTSKY database and the archive of the CoRoT 
Data Center). 

A significant number of stars with R < 16 is ex- 
pected to be giants and supergiants according to simu- 
lations (|Robin et al.l l2003l ) of the stellar p opulations in 
the fi el d of view. A s fully explained in [Peleuil et al.l 
ll2009h. ICabrera et al.l (|2009h . lAigrain et al.l (|2009l ) aiid 
ICarone et al.l (1201 lI V dwarfs and giants are separated in 
colou rs, using SED analysis techniques (see iDeleuil et al.l 
I2OO9I and references therein, for further information), in 
color-magnitude or colour-colour diagrams, combining the 
information of the EXO-DAT optical photometry with 
the near-infrared Two Micron All Sky Survey (2MASS) 



Point Source catalog jCutri et al.ll2003 ). Comparing the re- 
sults with simulations obtained with a model of synthetic 
stellar populat i ons of the Galaxy (the Besancon model, 
iRobin fc Crei^ 1 19861 : iRobin et al.1 bOOSl l. taking also into 
account the inhomogeneities of the galactic pl ane and the 
mean reddening of the fields (jRuphv et al.l[l997l l. they found 
that the two populations are quite well separated at the 
brighter end of the CoRoT magnitude range. For a field to- 
wards the center direction which is densely crowded, they 
estimated that less that 15% of the stars brighter than 
R = 14.5 are in the overlapping region between the two 
groups. There is some degeneracy for low-luminosity objects 
as the two populations mix and the separation becomes less 
reliable, particularly in the center direction. The dividing 
line between dwarfs and giants identified by Aigrain et al.l 
(200(|) in the R versus B — V colour diagram, as shown in 
Figure [31 was obtained from detailed checks with the pre- 
dictions of galactic models and constitutes an acceptable 
compromise, from a statistical point of view, between the 
loss of some (redder) bona fide dwarfs from the sample and 
the minimization of the residual giant contamination. To re- 
move most of the giants from our sample, we eliminated all 
the stars falling, in the R versus B ~ V colour diagram, to 
the right of this dividing line. 

Here we present the rotational analysis of the complete 
sample of 8341 white and chromatic (red, blue and green 
colours) hght curves selected with the procedure described 
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Figure 3. /{-magnitude versus B — V colour of the selected CoRoT targets in the anti-center {left) and center fields (right). Dots on the 
left of the continuous line correspond to stars identified as likely dwarfs, dots on the right of the line to likely giants. The cut between 
dwarfs and giants was empirically set as a straight line running from B — V = 0.7 and 1.1 at i? = 11 to _B — V = 1.2 and 1.2 at R =16 
for LRaOl and LRcOl, respectively. The larger fraction of giants in LRcOl is clearly visible. 



Table 1. Number of LRaOl (top) and LRcOl (bott om) ligh t curve s 
for main-sequence stars, selected following Aigrai rTet al.l l l2009t) , 
amongst the different colour indices and magnitude bins. 
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above (we eliminated 2934 LCs, classifying them as giants). 
For the aims of our work we did not take into account the 
colour information and we simply used the white hght curve 
obtained by summing up the fluxes in the three RGB bands. 
In particular we analyzed 5339 LRaOl hght curves, and 3002 
LRcOl light curves (Table [T]). 



We started from the N2 light curves in which most, 
but not all, instrumental effects have been corrected. The 
residual artifacts are: 

(i) long term drifts. These variations are partly due to 
pointing and to instrumental drift; 

(ii) residual periodic variations on the timescale of the 
satellite orbital period (1.7 h); 

(iii) discontinuities due to hot pixels, which often become 
unstable, causing further jumps through the run after the 
first jump (although the pipeline flags almost all of them 
and removes most). 

In order to mitigate the negative effects of these arti- 
facts for our analysis, we preprocessed the light curves in an 
iterative way. We implemented a routine which, as a first 
step, automatically analyzed all the LCs, detrending, rebin- 
ning to 2h and separating them in two groups according to 
the presence of discontinuities. As a second step, we visu- 
ally inspected the results and intervened case-by-case, when 
necessary, for further correction and selection. In particular: 

(i) To eliminate the long term drifts we detrended the 
light curves by fitting a third degree polynomial to the data 
(which works well for most of our LCs) and then dividing 
the original data points by this fitting. An example of the 
result can be seen in Figures [T] and (2] in which on the 
left we have the original light curve and the detrended one 
on the right. To ensure that no strange behaviour resulted 
from the process (in particular for the LCs with jumps) the 
results were inspected by eye, and eventually corrected for 
the residual instrumental variations with the appropriate 
fit decided on a case-by-case basis. In particular, for the 
LCs with jumps which we judged usable (cutting away the 
segments with discontinuities, for example) the detrending 
process was reiterated. 

(ii) We rebinned the data to two hours to smooth out 
the orbital period (1.7 h). Eliminating the signature of the 
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orbital periodicity reduces the effect of tlie outliers (mainly 
due to tfie satellite crossing of the South Atlantic anomaly). 

(iii) Many CoRoT LCs show sudden large flux discon- 
tinuities, and in a discrete number of cases each LC 
is affected by several "hot-pixel" events. If only one 
pixel is affected, this can be identified by comparing the 
flux in different color channels (when available). Differ- 
ent CoRoT use rs use different techniqu e s to deal with the 
jump problem jAl apini fc Aigrain' '2008'; 'Bord e et aLll2007l : 
ICarpano fc Fridlun d 2008; Ronnor ct al. 2008:), and there 
are several examples of employed strategies also for the 
"corr ection" of Kepler LCs with discontinuities l|Basri et al.l 
|2010| ). Nevertheless, since there is no simple way to cor- 
rect all the CoRoT LCs displaying jumps (due to the un- 
predictable gradual decay which follows each jump and the 
complicate behaviour of each LC affected by several jumps) 
without altering the astrophysical signal, we simply decided 
to retain only LCs with discontinuities smaller than a se- 
lected threshold, verifying by visual inspection that discon- 
tinuities do not affect detrending procedure. In order to rec- 
ognize the jumps we have first computed the distribution 
of the differences between fiuxes in adjacent bins and evalu- 
ated mean value and standard deviation, an- In Figure |4] we 
report a couple of examples. Since most of the LCs contain 
more than one significative jump, we utilize the 3 largest 
adjacent flux differences (which we call Mi , Mg and Ms, see 
Figure \^ , and define the parameters for identifying the two 
discontinuities as follows: CRj = (Mj — Mj+i)/(jd, with 
J=l,2. 

We selected all the light curves which have both CRi and 
CR2 lower than a selected threshold, which we set arbitrar- 
ily to 10, which means that we reject light curves with jumps 
that vary from the mean value by more than 10 standard 
deviations. There are several LCs, which did not fulfill our 
criteria for their selection while showing a clear periodicity 
by visual inspection, and we decided not to reject these LCs. 
Our selection criterion results in about 28% light curves re- 
jected for a total of 6001 LCs analyzed. The group of rejected 
LCs were inspected by eye, to evaluate if it was possible to 
use a shorter segment of the light curve, eliminating that in- 
cluding the discontinuity, to find a periodicity. We retrieved 
240 "cutted" LCs, which have been properly corrected for 
long period trend and have been analyzed serching a pe- 
riodicity, as well as the rest of the sample. In Figure |4] two 
examples of rejected light curves are shown. Top panels show 
a case with CRi > 10, and bottom panels a case with CR2 
> 10 and CRi < 10. 

Using lower thresholds, namely 5 and 3 standard devi- 
ations, we would reject about 47% light curves and about 
59% light curves, respectively, but we have verified that final 
results do not change. We conservatively decided to adopt 10 
a as threshold. Pathological cases with several small jumps 
and no periodicity are rejected by the autocorrelation anal- 
ysis (see Sec. 14. 2[) . 



4 ANALYSIS OF THE COROT LIGHT CURVES 

The analysis consists of two different methods: the Lomb- 
Scargle periodogram analysis, the autocorrelation analysis 
and a consistency check supported by the visual inspection 
of each light curve. Our final result is the compilation of a 



catalogue of periods for field main-sequence stars, towards 
the Galactic center and anti-center. From a detailed analysis 
we determine which of these periods can be attributed to 
rotation. 



4.1 Periodogram and folded light curves 

We searched the light curves for the presence of signif- 
icant periodicity using the L omb-Scargle Normalized Pe- 
riodogram (LNP) app roach jLombl 1 19761 ; IScarglel Il982l ; 
iHorne fc Bahunadll986l ). 

The periodogram provides an approximation to the 
power spectral density of a time series. The approach used 
makes periodogram analysis exactly equivalent to least- 
squares fitting of sine curves to the data. A function P{io) 
is computed from the data for a range of frequencies (which 
we have set from 0.01 to 10.0 d~^, to approximately match 
the minimum and the maximum searchable periods in the 
observational window) and the frequency that maximizes 
this function is considered as the most likely frequency in 
the data. Thus, we calculated the normalized power P{io) 
as function of angular frequency uj = 2tt v and identified 
the location of the highest peak in the periodogram. 

In order to decide the significance of the peak we have 
followed I Eaton et al.l (|l995l ). randomizing the temporal bins 
from the original light curve. By calculating the maximum 
power on a large number of randomized data sets, the con- 
version from power to False Alarm Probability (FAP) can 
be determined. In detail, we constructed 1000 light curves 
resampled from the original ones randomizing the position 
of blocks of adjacent t emporal bins (block length, 12 h) (e.g. 
iFlaccomio et al.ll2005l ). By shuffling the data we break any 
possible time correlation and periodicity of the light curve 
on time scales longer than the block duration. We calcu- 
lated the Scargle periodogram for all the randomized light 
curves and we compared the maximum from the real peri- 
odogram to the distribution obtained from the randomized 
light curves, at the same frequency, in order to establish the 
probability that values as high as the observed one are due to 
random fluctuations. Given some threshold FAP, we state 
that the detected candidate periodicity is statistically not 
significant if FAP > FAP,. The calculation we performed 
on CoRoT light curves led often to small FAPs, indicating 
that LNPs of our light curves present peaks that in most 
cases cannot be explained by pure stochastic noise, or non- 
periodic variability on time scales shorter than 12 h. In fact, 
if light curves present variations on time scale smaller than 
the size of the temporal block we used in the simulations, 
these variations will be still present in the simulated curves 
and will be not recognized as significant. We have chosen 
a bin size of 12 h as a reasonable compromise between the 
expected time scale of stochastic variations and the shortest 
expected periodic signal. 

In particular out of the total analyzed sample of 6241 
LCs, 4887 have significant periods (at FAP 1%; since we 
simulated 1000 light curves, the 1% FAP power is the power 
that was exceeded by the highest peak in 10 simulations). 
We show in Figure [S] some examples of the periodogram 
analysis. 
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Figure 4. Top left panel: Example of a light curve affected by a hot pixel event a few days after the start of the observations, causing a 
sudden rise in the measured flux followed by a gradual decay. Several "spikes" (t~ 2730 d, t~ 2760 d and t~ 2765 d) are present in the 
flux, due to high energy particle impact onto the CoRoT CCD. Top right panel: Distribution of the differences between adjacent fluxes: 
we indicated the mean value, the cr/j, the differences {D2,D3) between the most deviant values and the two rejection criteria Ci?j and 
CRs- The criterion used for rejection is CRi > 10. Bottom left panel: In this second example we underline the presence of two large 
jumps in the flux, at t^ 2705 d and t~ 2820. Bottom right panel: The criterion used for the rejection of this light curve is CRg > 10. 



4.2 Autocorrelation analysis 

Aliasing effects or residual effects due to the particular 
choice of temporal blocks used in the light curve's simula- 
tion may cause spurious per iods. To overcome th ese we used 
an autocorrelation analysis (|Box fc Jenkinsiri976l ). mainly to 
confirm the periodicity of the light curves, or to identify the 
correct period. Autocorrelation takes each point of the light 
curve measured at time t and compares the value of that 
point to another at time t + L. Points separated by L will 
be very similar if the data contained some variability with 
period L, thus the autocorrelation function will have peaks 
corresponding to periods of variability in the data. The au- 
tocorrelation r_L of a sample population X as a function of 
the lag L is: 

_ I]fro^~^(afc - x){xk+L - X) 

where x is the mean of the sample population X and N is 
the sample size, the quantity tl is called the autocorrelation 
coefficient at lag L. The correlogram for a time series is a 
plot of the autocorrelation coefficients as a function of L. 
A time series is random if it consists of a series of indepen- 
dent observations with the same distribution. In this case 
we would expect the to be statistically not significant for 
all values of L. We have chosen to adopt a 95% confidence 
level to select significant autocorrelation coefficient. Follow- 



ing this criterion a total of 3578 LCs present a significant 
autocorrelation . 

In Figure |S] we show two examples of autocorrelation 
plots for two stars of our sample. In this figure we can no- 
tice the particular behaviour of the top LC, the one with 
CoRoT ID 0102697871. The L-S periodogram found three 
significant periods (Pi = 47.5, P2 = 32.7, P3 = 24.9 d) at 
variance with the correlogram, in which no significant pe- 
riod is found, moreover the phase foldings performed with 
the three periods indicated by the L-S method show a very 
large scatter. We analyzed in detail this star in order to un- 
derstand its behaviour and discriminate between astrophys- 
ical and instrumental signals. In particular, we compared 
the evolution of the flux in the three different channels (red, 
green and blue) which were available for this star, and which 
show that instrumental signatures significantly perturb the 
LC. There is a data gap of about 4 days between 2939 and 
2943. A proton impact led to a reset of the Data Processing 
Unit (DPU) 1, which is responsible for data collection on 
the El CCD, on January 18, 2008 during SAA crossing. AU 
light curves in LRaOl originating from CCDl contain this 
data gap. The stepwise periodic changing between 2870 and 
2900 is not real but it is due to a peculiar combination of hot 
pixel of high and low amplitude affecting different channels. 
In particular, in T = 2858, 2867 two small amplitude hot 
pixels appeared in the blue channel, and two bigger events 
at T = 2870 and 2890, then two big ones affected the red 
channel at T = 2876 and 2882, the rest of the up and down 
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Figure 5. Examples of the detrended light curves (left) for two CoRoT sources (observed towards the Galactic center); periodograms 
(right) for the two LCs, the dotted curve superimposed on the LNPs is the threshold determined by simulations. Periods found with the 
periodogram analysis for the star with CoRoT ID 0100993180 (top) are not significant since they have Porig—LC < P simul—LC > while 
periods found for the star 0100889824 {bottom) are significant since they have Porig—LC > Psimul-LC- Note the very different vertical 
ajcis. 



behaviour is due to strong hot pixel events in only the blue 
channel. The red and the green channel display opposite 
trends (as usual) with the red channel affected by hot pixels 
in the first 30 days of the LC, thus also the large variation 
shown by the white curve seems to be instrumental. There 
are several light curves displaying a similar behaviour as the 
illustrated one, and we performed an analogous analysis us- 
ing colour channels, when available. In other cases, with no 
colour information, we had to trust the cross-correlation of 
the results obtained with L-S periodogram, phase folding, 
autocorrelation technique and visual inspection, but always 
following a conservative approach which minimizes the num- 
ber of spurious periods in the final catalogue. 

4.3 Consistency check and visual inspection 

After having derived periods both from the periodogram 
and the autocorrelation analysis, we have checked the con- 
sistency of the results obtained with the two methods. In 
none of the cases in which the Lomb-Scargle method does 
not find any significant period, the autocorrelation method 
finds one, that implies that the latter is more conservative 
in identifying periodicity. A total of 1457 stars have consis- 
tent periods according to both methods. In 1152 cases we 
have periods significant for the Lomb-Scargle analysis but 
not significant for autocorrelation. Since we want to derive 
a sample with confirmed periods, we have discarded these 
light curves. In 2121 cases both methods give discrepant sig- 
nificant periods. We have examined by eye the phase folding 



with both periods, choosing as final period the one that re- 
produce a smoother behaviour. In many of these cases the 
rejected period is an alias, and the folding helps in discrim- 
inating the original period (Figure [71 left). This approach 
has made possible to determine periods for 521 LCs , in the 
remaining cases the folding gave unsatisfying results (very 
scattered plots) and therefore we discard them from the pe- 
riodic LCs. We want to stress that there may be a complex 
mixing of several phenomena (i.e. chaotic activity, fast evo- 
lution of spots, multiple periods which interfere with one an- 
other, residual instrumental effects) which makes the identi- 
fication of a clear periodicity extremely difficult, on the basis 
of the photometric behaviour only. For these stars it would 
be useful to analyze both LCs collected at different times 
and spectroscopic observations. In a large fraction of these 
discarded periods, we noticed the presence of several small 
jumps, that are not identified by our procedure described 
in Sec. |31 that strongly affect the analysis producing spuri- 
ous periods. An examples of such cases (spurious period) is 
shown in Figure [7] (right). 

In total the described procedure has allowed us to de- 
termine periods for 1978 stars (1457-1-521). 



5 RESULTS: IDENTIFICATION AND 

ANALYSIS OF ROTATIONAL PERIODS 

In order to study the rotational properties of our sample we 
need to identify the eclipsing systems and pulsators. Eclips- 
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Figure 6. Two examples of light curves {left panels) and relative correlograms (right). In the autocorrelation plots we marked the 95% 
confidence interval (dotted horizontal lines) and the first (solid vertical line) and second period (dashed vertical line) found by the 
periodogram analysis (the values of the first two periods yielded by the L-S periodogram are indicated in the top right of each plot). 



ing stars are easily identified by eye from the narrow eclipsed 
region in the folded curve (see an example in Figure |8l left 
panels). We identify a total of 82 eclipsing systems. How- 
ever in a fraction of them it is possible to derive both the 
orbital and rotational periods. In these cases we report both 
periods in Table 2. 

A further step is the identification of pulsators (see for 
example Figure [H] right panels). 

The distinction between pulsations and rotational 
modulation is difficult to make solely based on the fight 
curve. To distinguish between variations due to pulsation 
and rotation caused by spots we adopted a method which 
exploits both the stability of the detected periods for long 
time and the small flux d ispersion in the folded curve 
l|Schmidt et al.ll2004 lioOSal lbl). 

Various observable parameters have been proposed 
to distinguish be t ween different types of pulsators. 
IWeldraJte fc BavlissI (|2008l ) define pulsators as stars dis- 
playing regular pulsation brightness variations with periods 
less than 2 days. The majority of these pulsators displays 
regular radial sinusoidal variations. Longer period pulsators 
are defined as stars displaying regular sinusoidal or slightly 
sawtooth variations with various periodicities greater than 
an arbitrary chosen lower limit of 2 days. These stars are 
likely giant stars of late spectral types. Apart from the 
definition used and the particular form of the LCs, one 
distinctive characteristic of pulsators is the longer ter m 
stability of the pulsation (|Schmidt et all |2004 l2005al lbl: 
iDebosscher et al.ll2009. ). 

We adopted a criterion based on the hypothesis that 



pulsator LCs are much stabler than rotation LCs. Then we 
divided the phase interval in 40 bins and we calculated the 
median flux value for each of these bins. Through linear 
interpolation of these median values we derived the median 
flux as a function of phase (see Figure |9} and we calculated 
for each phase value the ratio between the difference of the 
instantaneous flux and the median phase folded lightcurve 
at the corrisponding phase and the peak-to-peak amplitude 
variation of the LC, Fluxi-Fluxi„terpi/(Fluxmax-Fluxmi„). 
We identify a star as a pulsator if more than 80% of the 
Fluxi-Fluxi„ie,pj/(Fluxmaa:-Fluxmin) ratios are lower than 
0.15. These thresholds have been chosen after a detailed 
inspection of all the LCs, identifled with pulsators by 
eye (LCs very regular with very small dispersion), and of 
a comparison sample of clearly non-pulsating LCs. The 
described method retrieves about 90% of the pulsators 
identified by eye. To be conservative, we add the remaining 
10% LC to the pulsator sample. We have decided to 
not increase the chosen threshold to include them, since 
otherwise we would include a significant number of LCs 
with evident characteristics of rotators (more irregular and 
higher dispersion LCs). A total of 169 pulsators has been 
fou nd. 111 of which h a ve per iods smaller than 2d. 

IDebosscher et al.1 (120091) performed an extensive au- 
tomated classification of variable stars for the first four 
CoRoT exoplanet fields, stating that we do not expect to 
find many classical radial pulsators while more candidate 
multiperiodic, non-radial, pulsators are found, due to 
the selection of the CoRoT exoplanet observation fields 
biased towards cool main-sequence stars. Several stars 
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Figure 7. In the left panels is shown an example of a LC {top) for which the most significant period derived with the Lomb-Scargle 
method P(LS) is an aUas of the period indicated by the autocorrelation analysis P(AC), as can be seen in the relative phase foldings 
{middle and bottom panels, respectively). In the right panels is shown an example of LC {top) affected by several jumps (smaller than 
the 10 (T threshold). In this case the two methods find discrepant periods. On the basis of the analysis of the relative phase foldings 
{middle and bottom panels) we decided to discard these periods. 




Figure 8. Two examples of light curves and folded LCs (first Lomb-Scargle period) for an eclipse star (ID: 0102588918, left panels) and 
for a pulsator (ID: 0102592190, right panels) . 



identified as pulsators with our method are probably 
giant stars, contaminating our sample. Amongst selected 
pulsators we recognized several double-mode RR-Lyrae, 
few low-amplitude Cepheid pulsators and a few RR-Lyrae 
displaying a clear Blazhko efTect (the LC varies in height 
and shape from cycle to cycle due to still poorly understood 
processes which afTect the regular pattern). 

We compared our cla ssification of pulsators with 
the one made by iDebosschc r ct al.. (|2009t ) for the com- 
mon objects and we found that, using our method, we 
always succeed in making the same classification as the 



iDebosscher et al.l (|2009l ) one. 

Amongst the LRaOl identified pulsators we found 16 
object which were also observed during the IRaOl (the 
initial run toward the galactic anti-center, with a data 
lenght of 55 days, which partially overlaps with the LRaOl 
field). As a further verification, we analyzed the LCs for 
these objects, the periods found for the two sets of data 
(IRaOl and LRaOl) are the same, and we classify these 
16 IRaOl objects also as pulsators, following our method. 
There is a gap of « 206 days between the two observing 
runs, thus we phased the two sets and we found that the 
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Figure 9. Interpolation curve of the median folded flux values 
for a candidate pulsator. 



pulsations are stable during this longer range of time. 
Such long time stability of peri ods is highly imp r obable 
for spot-like modulations. As iDebosscher et aP (|2009l ) 
admitted, rotation can produce LCs that are difficult to 
distingush, for example, from the typical skew-symmetric 
Cepheid LCs, nevertheless we need spectral information 
and LCs with a longer total time span to verify whether 
the detected periods remain stable over time. 

We remain with 1727 stars whose modulation is inter- 
preted as due to rotation. All the derived periods and their 
classification are reported in Table (2] The table shows a 
portion of the catalogue: the first and second columns re- 
port the CoRoT and 2MASS ID of the stars; the third and 
fourth their right ascension and declination; the 5th column 
reports the derived periods; the 6th and 7th are flags indicat- 
ing the method(s) adopted to derive periods and the nature 
of the stars (pulsator, eclipse star, etc), respectively and the 
8th is a link to the plots of the detrended light curve and of 
the folding. The 9th column reports the rotational periods 
of eclipsing binaries, when derived. 

The distributions of rotational periods for F, G and K- 
type stars (both Galactic anti-center and center) are shown 
in Figure 1101 the distribution of periods of F stars in the 
center direction is missing because of their paucity, only 4 
stars and none with significant period. 

The distributions shown in Figure [TOl are bimodal, with 
two peaks probably associated to two distinct populations, 
one (presumably young stars and binaries) peaking at very 
short rotational periods (Prot < 5 - 10 days), and the other 
(older) peaking at longer rotational periods, with a gap be- 
tween about 15 and 35 days. 

The presence of a large number of stars with short ro- 
tational periods is fully compatible with the presence in the 
solar neghborhood of a sample of young stars, both towards 
the center and the anti-center of the Galaxy. Stars in the 
short period peak in Figure [10] have periods shorter than 
typical for old single stars. They are likely dominated by 
young stars that have not had the time to spin down with 
perhaps a fraction of older tidally locked binary stars. In 
fact, by visual inspection of the light curves we have been 
able to detect some binary systems of different type (eclips- 
ing, contact and detached binaries). 



0.30 



0.25 



S 0.20 



F,G-type 



O. 

o 

a 

3 



0.15 



0.10 



0.05 



0.00 




0.20 



20 40 60 80 100 

P_rot (Days) 

K-type 



0.15 



c 

(9 

O. 

o 

a 

3 



0.10 



0.05 



0.00 



K.,(254) 
K (396) 



11 



10 15 20 25 



20 



40 60 
P_rot (Days) 



80 



100 



Figure 10. Normalized distributions of rotational periods for the 
sample analyzed of F and G-type stars {top panel) and K-type 
{bottom panel) stars towards the Galactic anti-center and center. 
The distribution of periods of F stars in the center direction is 
missing because of their paucity (4 stars). The numbers reported 
indicate the number of light curves for each type. 



6 DISCUSSION 

As described in the introduction, rotational period may 
be used as a proxy for age for single dwarfs of a given 
mass. In particular this is possible for stars older than 10^ 
years when solar type stars have converged to a single value 
rotation-age relation. This has been confirmed by a num- 
ber of studies of rotational rate in intermediate-age open 
clusters (e.g. the 150 Myr old M35. iMeibom et alj [ioool . 
the 250 Myr old M34, llrwin et al.l |2006| . the 550 Myr old 
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Table 2. Catalogue of periods for CoRoT main-sequence stars towards the Galactic center and anti-center (the complete table is available 
online). 
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Figure 11. Rotation period vs B — V derived from CoRoT light 
curves for solar-type sta r s com pared to the relations derived in 
iMamaiek fc Hillenbrandl ||2008| ') for 130 Myr (Pleiades) and for 
625 Myr (Hyades). 



M37. iMessina et a l.' '2008*. the 600 Myr old Coma cluster, 
Cameron et al..,200& and the w ell studied 625 Myr 
old Hyades cluster, Radick et al.l Il987l ) which have shown 
that the rotation of F and G stars at the age of 150 Myr 
is already single-val ued and only depend ent on color (that 
is, mass) and age, (jMeibom et al.ll2009h . while for K-type 
stars this happens at the age of the Hyades, at which rapid 
rotators are still present amo ng M-type stars. W e refer the 
reader to the recent work of llrwin et al.l (|201ll ) (and ref- 
erences therein), for a detailed discussion on rotation and 
angular momentum evolution of M-type field main-sequence 
stars from the MEarth transit survey. They measured a wide 
range of rotation periods for 41 M dwarfs (from 0.28 to 154 
days), finding that kinematically young (thin disk) objects 
rota te faster than the kinematically old (thick disk) objects. 
iBarned |20o3) used literature data to derive the gy- 



rochronology relation. According to him, the rotation pe- 
riods evolve with age as: 

P{B~V,t) = f{B-V)g{t), 

showing a dependence on both age and color. 

Since our sample includes very few dM stars, we may 
safely consider that stars slower than the Hyades rotation- 
mass relation, are older than 600 Myr. Faster stars are cer- 
tainly younger. Therefore we used the Hyades rotation-mass 
relation (Barnes 2007; Mamajck & Hillenbrand 2008), as a 
dividing line between stars younger and older than 600 Myr, 
as shown in Figure 1111 Gyrochronology relations for aged 
main-sequence stars rely only on the Sun and few other 
stars, thus, even if they are older, we decided to exclude 
rotational periods greater than 30 days from the plot. 

As can be seen in the Prot vs B — V plot, more than 
500 stars rotate faster than the mean rotation rate expected 
for the Hyades, and we consider these stars as dominated 
by a young population (^ 600 Myr). Our sample suffers of 
several biases preventing us to perform a comparison with 
the observations and the predictions from Galactic models. 
In fact the original sample was selected to maximize the 
probability to detect planetary transits, on top of that we 
discarded from our analysis a fraction of targets affected by 
instrumental artifacts, and finally we have not been able to 
identify rotational periods in a significant fraction of stars. 
There are several reasons for that, the most obvious is that 
we need a minimum level of activity to detect a minimum ro- 
tational modulation. This will bias the final sample against 
the more quiet stars, likely the oldest ones. On the other 
extreme, chaotic activity, characteristics of very young stars 
may prevent the identification of a clear periodicity. 

Stellar rotation is one of the best proxy for young 
main sequence stars, since their photometric and spectro- 
scopic properties do not differ significantly from those of 
the oldest field stars. However this population may be con- 
taminated by close binary systems, which have maintained 
higher rotation rates due to tidal synchronization. We can 
estimate the fraction of binary stars in our s amples from 
the w ork of iDuguennov fc Mavorl l|l99lh and iMavor et al.l 
(| 19921 ), on the statistics of binarie s in the solar neighbour- 
hood. IDuguennov fc Mavorl l|l99ll ) found that the fraction 
of binaries among G-type stars in the solar neighbourhood 
is about 60%, 6% of which is constituted of short period 
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binari es, those of interest here (P < 20 days') . [Mayor et al] 
found that the analogous fraction in K-type stars 
is about 45%, with a contribution of about 3% from short 
period binaries. These estimates imply that binaries may ac- 
count only for a small fraction of our fast rotator sample, 
then we can safely identify the bulk of our fast rotators as 
composed mainly of single stars. 

Previous studies based on the stellar content of shal- 
low and intermediate X-ray surveys have revealed a dis- 
crepancy between model star counts based on constant 
rates of star formation and observed star counts for ac- 
tive stars. In particular, an excess of X-ray bright FGK 
stars has been observed (iFayata et al. l ll988l : ISciortirio et al.l 
1 19951 : iLopez- Santiago et al.l 20071 ). for example, by the 
Einstein Extend ed Medium Sensitivity Survey (EMSS, 
iGioia et all Il990l), the XMM-Newton B right Serendipitous 



Survey fXBSS. iDella Ceca et al.] |2004) and the ROSAT 
North Ecliptic Pole survey (iHenry et al.ll200ll '). These late- 
type stars appear to be young, as shown by their lithium 
abundance (Favata ct al. 1 993t) , although a fraction may be 
active binary systems that h ave a yellow dwarf as a primary 
jLopez-Santiago et al.l [20071 V Data are consistent with the 
presence of a recent event of star formation in the last 10* 
years, while in the range 10^-10^" years, the assumption of 
a model with a c onstant rate of st ar formation agrees with 
the observations ([Affer et al.ll2008l '). 

Analogously several young (8 - 50 Myr) moving groups 
have been identified both i n the northern and southern sky 
llZuckerman fc Sond [2004'. 'Lop ez-Santiago et all [ioo^ and 
reference therein), suggesting that the history of the star 
formation close to the Sun may be very complex with a few 
recent bursts. 

In light of the mentioned suggestions, reported in the 
literature, of a recent event of star formation in the solar 
neighbourhood, it is plausible that the fast rotators identi- 
fied with our analysis trace this young population. Spectro- 
scopic observations of these stars would clarify their nature 
and provide a further test of the presence in the solar neigh- 
bourhood of the young population indicated by the activity 
surveys. 



der to maximize the probability to detect exoplanetary sys- 
tems eliminating stars affected by crowding. The biases in- 
troduced by this choice are not obvious. Furthermore the 
data analysis procedure has discarded a non- negligible frac- 
tion of light curves due to instrumental artifacts. Finally, 
as discussed in the previous section, we may derive periods 
only in stars with clean light curves and with amplitude 
variations measurable in our observations. 

However our analysis has allowed us to identify pul- 
sators, binaries and, among a huge number of field stars, 
a sample dominated by young stars with small contamina- 
tion by binary systems. The nature of our selected sample of 
young stars is consistent with the nature of the population 
identified by X-ray limited surveys that shown the presence 
of an excess of young stars in the solar neighborhood. 

Future spectroscopic follow-up observations are needed. 
This will also allow us to derive chemical abundances, in 
particular lithium, and radial velocities, as well as to obtain 
information from several activity indicators, for instance the 
emission of chromospheric features such as Ca II H and K. 

Spectroscopic observations will allow us to obtain an 
independent estimate of the rotation of our stars through 
the measure of the projected rotational velocity, «sini. This 
measure may be extended to a sample of stars for which we 
were unable to measure the period in particular for stars 
with chaotic LCs, in order to better understand the biases 
in our period determinations (i.e. to show if these stars are 
mainly slow rotators). 
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7 SUMMARY 

We have analyzed the light curves observed by CoRoT in two 
fields at different direction in the Galaxy, in order to derive 
stellar rotational periods. Our original sample of 8341 light 
curves has been reduced to 6241 due to instrumental arti- 
facts. Light curves have been analyzed both with the Lomb- 
Scargle periodogram and with the autocorrelation methods. 
After a careful check of the consistency among the results 
and a visual inspection of all the original and folded light 
curves we retain 1978 periods. After having identified eclips- 
ing systems and pulsators we identify 1727 rotational peri- 
ods. 

Taking advantage of the more recent version of the gy- 
rochronology relations, we have been able to identify more 
than 500 stars in our sample which are likely younger than 
600 Myr. We cannot use our rotational period distribution 
to infer the age distribution in the solar neighborhood since 
we cannot assess the completeness of our sample. In fact 
the original sample observed by CoRoT was selected in or- 
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